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Abstract 

We discuss the reconstruction of neutrino flavor neutrino at a distant source in the very high en- 
ergy regime. This reconstruction procedure is relevant to the confirmation of detecting cosmogenic 
neutrinos, for example. To facilitate such a reconstruction, it is imperative to achieve effective 
flavor discriminations in terrestrial neutrino telescopes. We note that, for energies beyond few tens 
of PeV, a tau-lepton behaves like a track similar to a muon. Hence, while it is rather challenging to 
separate from in this case, one can expect to isolate from the rest by a distinctive shower 
signature. We present the result of flavor ratio reconstruction given the anticipated accuracies of 
flavor measurement in neutrino telescopes and current uncertainties of neutrino mixing parame- 
ters. It is shown that the further separation between and Vr events does not improve the flavor 
reconstruction due to the approximate — v^- symmetry. 
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I. INTRODUCTION 



High energy neutrinos with different flavor ratios can be produced from various astro- 
physical environments or production processes. For instance, a typical hadronic interaction 
produces neutrinos through tt^ ^ /i^ — )■ decay chains and generates a neutrino flavor 
ratio, (poi^e) '■ (poi^fi) '■ 0o(^t) = 1:2:0 [1]. This type of source shall be referred to as pion 
source hereafter. When the source is hidden behind a strong field [2] or a dense matter jsl, 
the muon produced by the pion decay may lose a significant amount of its energy before 
its decay j4|. This type of source is referred to as muon-damped source since the neutrinos 
produced by the decay of secondary muon give negligible contribution to the total neutrino 
flux. Hence the resulting neutrino flavor ratio is (poi^e) '■ 4'o{i^fi) '■ 0o(^t) = 0:1:0. 

The flavor ratio of neutrinos are altered by oscillation effects while neutrinos propagate 
from the source to the detector on the Earth. Since the distance from an astrophysical source 
to the Earth is generally much longer than the neutrino oscillation length, the observed 
neutrino flavor ratio is fully determined by three neutrino mixing angles and the CP phase. 
In this work, we aim to reconstruct the neutrino flavor ratio at the source from the observed 
flavor ratio on the Earth. As we shall see later, the method for such a reconstruction depends 
on the neutrino energies. We shall focus our discussion on the highest energy regime, which 
is an extension to our earlier work 



II. CLASSIFICATION OF NEUTRINO EVENTS AT DIFFERENT ENERGIES 

Neutrinos interact with matters to produce observable signals. The major channel for 
such interactions is the charged- current (CC) interaction, ui + N — > I + X, where / is the 
lepton associated with ui and X denotes the hadronic states. The sub-dominant channel is 
the neutral-current interaction (NC), vi + N — y vi + X. In Fig. [Hand Table HI we summarize 
different types of neutrino induced events and their detectable energy ranges. 

Type-A event in Fig. [1] is an electron production through v^. CC interaction. The electron 
has a large interaction cross section with the medium and produces a shower within a short 
distance from its production point. Type-B event is a muon produced by CC interaction. 
Contrary to the electron, a muon can travel a long distance in the medium before it loses all 
its energy or decays. The muon range in ice is more than 10 km for E^y = 1 PeV (10^^ eV). 
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particle major processes signal type symbol in Fig(T] 

e EM shower shower A 

energy loss track B 

t{Ei, < 3.3 PeV) CC int. and r-decay shower C 

r(3.3 PeV< E^, < 33 PeV) CC int. and r-decay 2 separate showers D (double-bang event) 

t{Ei, > 3.3 PeV) energy loss and decay track and shower E (lollipop event) 

t{Ei, > 3.3 PeV) CC int. and energy loss shower and track F (inverted lollipop event) 

r(£'i^ > 33 PeV) energy loss track G 

X hadron shower shower H 

TABLE I: Different types of neutrino induced events. 

Hence, above this energy, there is hardly any decay of muon occurring within the fiducial 
volume of the detector, which is about a few km^. A muon does, however, lose a small 
fraction of its energy and emits dim lights so that only those optical detectors which are 
near to the muon track can be triggered. As a result, a muon produces a track-like signal. 

The i^T-induced events are listed as types C-G where the tau lepton produced by Ur CC 
interaction behaves differently at different energies for a fixed detector design. For a neutrino 
telescope such as IceCube [a], the distance between each string of optical detectors is 125 
m, which corresponds to the decay length of a 2.5 PeV tau lepton. Such a tau lepton could 
be produced by the CC interaction of a Ur with E,y = 3.3 PeV. Therefore, for a Ur with 
an energy significantly smaller than this, the separation between the first hadronic shower 
produced by CC interaction and the second shower produced by the tau-lepton decay is 
too small to be resolved. Such an event is classified as type C. For > 3.3 PeV, one can 
resolve the above double-bang event (classified as type D) until the separation of two showers 
exceeds the effective size of the detector. Such a size is estimated to be the sum of IceCube 
dimension (f« 1 km) and two extinction lengths of optical photons in ice (r^ 250 m), which 
corresponds to the decay length of a tau lepton with Er = 25 PeV. The average energy 
of Ur capable of producing such a tau lepton is around 33 PeV. Hence the configuration of 
IceCube detector determines the observable energy range for the double bang event to be 



3.3 PeV < E„ < 33 PeV 



For an under-sea experiment, such as KM3Net 10|], the 



observable energy range for the double bang event is similar. 
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Type-E event is referred to as the lollipop event. In such an event, a high energy tau 
lepton enters the detector and decays within it, producing a track signal followed by a shower. 
The probability for observing a lollipop event increases with the neutrino energy, and it is 
about 5 X 10~^ for Ey = 1 EeV fll|. Type-F event is the inverted lollipop which consists of 
a hadronic shower from Vj. CC interaction and a subsequent tau-lepton track. Both muons 
and tau leptons produce inverted lollipop events and it is not easy to separate them. Type- 
G event is a through-going tau-lepton track which is produced by CC interaction with 
Ejj > 33 PeV. Finally type-H event is induced by the neutral-current interaction. 



III. FLAVOR DISCRIMINATIONS FOR DIFFERENT NEUTRINO ENERGIES 

Although our focused energy range is at Ei, > 33 PeV, it is helpful to review the flavor 
discrimination in the lower energy ranges. 



Flavor discrimination for Ei, < 3.3 PeV 



In this energy range, type-C events can not be separated from type-A and type-H events 
since the two showers in type-C events can not be resolved. Hence one can only distinguish 
muon track event (type B) from shower events (type-A, C and H). In IceCube, such a 



distinction can be done effectively 



which is useful for deducing the flux ratio 



0(z/^)/(0(z/e) + 0(z/t-)) [12|. We note that 4>{i^^) in the numerator contributes to both type-B 
and type-H events. On the other hand, 0(t'e) and 0(^'t) contribute equally to type-H events if 
4>{'^e) = My-r)- Furthermore, type-A and type-C events occur with the same scattering cross 



section 



13| . This explains the flux combination 0(i^e) + 0('^r) appearing in the denominator. 



B. Flavor discrimination for 3.3 PeV < £"1^ < 33 PeV 

In this energy range, one can detect the type-D and type-E v-r events (double bang and 
lolhpop). Hence it is also possible to deduce the flux ratio = (j){ve) I 4>{i-'t) [l^ in addition 



to B}. However, the double bang and lo^ 



with is large. In an earlier paper js], we demonstrated that a large number of events is 
necessary for lowering down the errors of E} and to the point that one can distinguish 



lipop events are both rare so that the error associated 
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FIG. 1: Different types of neutrino-induced events. Dashed lines and solid lines correspond to 
paths of neutrinos and leptons respectively. The ellipsoids are showers. The energy range suitable 
for detecting each type of event is listed in Table HI 



Condition 1 : < 33 PeV Condition 11 : > 33 PeV 

TABLE II: Flavor discrimination variables for different energy ranges 
the pion source from the muon-damped source. 



C. Flavor discrimination for Ey > 33 PeV 

In this high energy regime, the tau-lepton range becomes long enough so that a tau 
lepton could pass through the detector fiducial volume without decaying. In this case, 
the tau-lepton loses its energy just like a muon does and the signal appears like a track 
event [l^. Thus, from an experimental point of view, one should classify such a signature 
as a track event (type G). In this energy range, there are also type-E and type-F events 
where tau leptons also behave like tracks. It is clear that the discrimination between tau- 
lepton and muon tracks are challenging, i.e, it is non-trivial to measure the parameter S^^ = 
0(z/^)/0(z/t-). However, one can expect to measure the fraction R^^ = 0(z/e)/(0(t'^)+0(z/T-)) 
by identifying the shower signature from the electron. In fact, can also be separated from 



other flavors in neutrino telescopes based upon radio detection technique la, llTI] due to 



Landau-Pomeranchuk-Midal effect 
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2Q|. 



As a short summary, we present in Table HTl the appropriate flavor discrimination variables 
for Ey < 33 PeV and > 33 PeV respectively. 



IV. THE RECONSTRUCTION OF SOURCE FLAVOR RATIO FOR > 33 PEV 



The evolution of neutrino flavor ratio from the source to the Earth is given by the prob- 
abihty matrix P^,^ such that 

(Pi^^a) = ^PapM'^l3)- (1) 

Each matrix element Pa/B is a function of neutrino mixing angles and CP violation phase, 
(pl^a) is the flux of neutrinos of flavor a on the Earth, and (poii'/s) is the flux of neutrinos 
of flavor P at the source. We assume the neutrino propagation distance is sufficiently large 
so that each P^/s depends neither on the neutrino mass-squared differences Am'fj nor on the 
neutrino energy E,^. 



A. The flavor reconstruction with only R^^ measured 

If one only measures R^^, the reconstructed range for neutrino flavor ratio at the source 



is given by the following fltting formula 21 1 



X'^=I%-^V+ E xL, (2) 



ij=12, 13,23 



where (T^ii^^ = (AP"/P^Pg^p. The functions Xg^^ taken from Ref. [22| while 



xi is taken from Ref. 23|. In the above references, the functional dependence of each Xe . 



on sin^ 9ij is given. The suffix "th" indicates the theoretical predicted values which depend 
on the source neutrino flavor ratio and the neutrino mixing angles. The suffix "exp" indicates 
the experimentally measured values which are generated by the true neutrino flavor ratio 
at the source and the best-flt values of neutrino mixing angles. The best-flt values and the 
allowed la and 3a ranges of mixing angles are given by 

„;„2/i _ ^ on/1+0-022,0.066 „\„2 n _ ^ c;+0-07,0.17 „;„2 a _ n m +0-009 /q\ 

sm Ui2 - U.dU4_o 054, sm 6^23 - U.5_o.o6,o.i4> sm 6^13 - U.Ul_o,oo65 l-^J 

where the la range for ^13 used in our analysis is that associated with normal hierarchy and 
the 3a range for ^13 in the same mass hierarchy is sin^ Oi^ < 0.035 [24,]. We take the CP phase 
6 = n for generating Pg^p- We have found that other 6 values do not produce noticeably 
different results. We consider all possible neutrino flavor ratios at the source for calculating 
X^- Since we have taken P^^p as those generated by input true values of initial neutrino 
flavor ratios and neutrino mixing parameters, we have Xmin — occurring at these input 
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FIG. 2: Reconstructed flavor ratios for an input pion source with AR^^/R^^ = 15%. The dark and 
light shaded areas denote la and 3cr ranges respectively. 

true values of parameters. The boundaries for la and 3a ranges of initial neutrino flavor 
ratios are given by Ax^ = 2.3 and Ax^ = 11.8 respectively, where Ax^ = ~ ix 



2^ ^ 2 

min A 



in this analysis. 

Let us first take the pion source as the input true source and consider its reconstruc- 
tion. The accuracy for measuring i?" is taken to be Ai?"/i?" = 15%. The possibility for 
measuring R^^ in such an accuracy will be discussed in the next section. The result for the 
reconstruction of neutrino flavor ratio is shown in Fig. [2l It is seen that the muon-damped 
source can be ruled out at the la but not on the 3a level. This result shows how well one 



can discriminate between cosmogenic neutrino flux (a pion source) 27 



interactions 
source 



31 



29 



arising from GZK 

30| and the high energy tail of astrophysical neutrino flux (a muon-damped 



32l |) when the former type of neutrino flux is detected. 
We next take the muon-damped source as the input true source. Such a neutrino flavor 
ratio occurs at the high energy tail of astrophysical neutrino flux as just mentioned. Once 
more we take AR^^/R^^ = 15%. It is seen from Fig. [3] that the pion source can be ruled out 
at the la but not at the 3a level from the flavor reconstruction. 
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FIG. 3: Reconstructed flavor ratios for an input muon-damped source with Ai?"/-R" = 15%. The 
dark and hght shaded areas denote la and 3a ranges respectively. 

B. The flavor reconstruction with both R^^ and 5^^ measured 

Although it is very challenging to measure S^^, it is of interests to see whether or not such 
an effort is useful for improving the reconstruction of neutrino flavor ratio at the source. Let 
us take a very optimistic scenario that the error of S^^ is comparable to that of R^^. The 
statistical analysis can be performed by adding S^^ contribution to Eq. ([2]). 

The reconstructed neutrino flavor ratio for an input pion source is presented in Fig. H] 
where AR^^/R^^ = 15% and AS^^/S^^ related to the former by Poisson statistics. This 
result is in fact comparable to that in Fig. |2] where only R^^ is measured, i.e., the further 
measurement of S^^ does not tighten the constraint on the initial neutrino flavor ratio. 



This can be understood by the approximate z/^ — Ur symmetry 33|, |34] which makes 5*^^ 
always rather close to unity no matter what the initial neutrino flavor ratio is. Hence the 
measurement of S^^ does not help to constrain the initial neutrino flavor ratio. We can 
demonstrate this further by taking the muon-damped source as the input true source. With 
AR^^/R^^ = 15% and AS^^/S^^ related to the former by Poisson statistics, the result for the 
flavor reconstruction is presented in Fig. |3 Comparing with Fig. [31 one can also see that 
the measurement of 5*^^ does not help to constrain the initial neutrino flavor ratio. 
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FIG. 4: Reconstructed flavor ratios for an input pion source with Aii"/i?" = 15% and AS^/S"" 
related to the former by Poisson statistics. The dark and light shaded areas denote the la and 3a 
ranges respectively. The result shown in this figure is comparable to that in Fig. [2] where only R^^ 
is measured. 

C. Comparison to the flavor reconstruction for E,^ < 33 PeV 



In this energy range, the appropriate flavor discrimination variables are and respec- 
tively as summarized in Table II. It has been shown in Ref. jsj that measuring significantly 
improves the constraint on the initial neutrino flavor ratio obtained by measuring alone. 
Taking pion source as the input true source, the effect of measuring is presented in Fig. Ei 
The left panel of the figure shows the reconstruction result obtained by measuring only R^ to 
the accuracy of AR^/R^ = 15% 35|]. At 3a level, any initial neutrino flavor ratio is allowed. 
The right panel is the result obtained by measuring both R^ and 5*^ with AR^/R^ = 15% 
and AS^/S^ related to the former by Poisson statistics. Clearly the further measurement of 

drastically improves the constraint on initial neutrino flavor ratio. 
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FIG. 5: Reconstructed flavor ratios for an input muon-damped source with AffVi?" = 15% and 
AS^^/S^^ related to the former by Poisson statistics. The dark and light shaded areas denote la 
and 3(7 ranges respectively. The result shown in this figure is comparable to that in Fig. [3] where 
only R^^ is measured. 




FIG. 6: Reconstructed flavor ratios for an input pion source in the energy regime Ei, < 33 PeV. 
The dark and light shaded areas denote the la and 3a ranges respectively. The left panel is the 
result obtained by measuring only to the accuracy of AR^/R^ = 15%. The right panel is the 
result obtained by measuring both R^ and with AR^/R^ = 15% and AS^/S^ related to the 
former by Poisson statistics. 
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FIG. 7 


: Reconstructed 


flavor ratios for an input pion 


source witfi Ai?"/i?" = 15% and 



uncertainties of neutrino mixing angles reduced to 50% of those listed in Eq. ([3|). The dark 
and light shaded areas denote la and 3cr ranges respectively. The result shown in this figure is 
comparable to that in Fig. [2] 

V. DISCUSSIONS AND CONCLUSIONS 



To reconstruct the source flavor ratio of ultrahigh energy neutrinos, it is crucial to have 
sufficient event numbers and an effective way of separating i/p, from and u-r- For > 33 



'eVjOne expects that the cosmogenic neutrino flux 



27 



GRB 



arising from GZK interactions 



29l . l30 [should dominate over those fluxes originated from astrophysical sources, such as 



36 



371]. It has been shown that the newly proposed Askaryan Radio Array [17] can 



detect roughly 50 cosnaogenic neutrino events in 3 years for baseline flux models such as 



those discussed in Ref. [28j. This event number implies a nearly 30% accuracy (statistically) 
in 3 years or a 15% accuracy in a decade of data-taking for determining the flux ratio R^^. 
Certainly the efficiency of flavor discrimination also affects the accuracy for determining R^^. 
Further studies are needed on this aspect. 

It is of interest to investigate whether or not a better knowledge of neutrino mixing 
parameters will improve the constraint on source flavor ratio. Taking the pion source as 
the input true source and considering the uncertainties of neutrino mixing parameters being 
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half of those hsted in Eq. ([3]), we obtain the reconstructed flavor ratios as shown in Fig. [71 
One cannot see noticeable difference between this result and that of Fig. |2l It is clear that 
the constraint on source flavor ratio solely depends on the accuracy of measuring R^^. 

In summary, we have argued that the flux ratio R^^ = (f){Vf.) / + 0(z^t-)) is a suitable 
variable for neutrino flavor discrimination for > 33 PeV in water (ice) Cherenkov and 
radio wave detectors. In view of recent development in radio wave array, we studied the 
reconstruction of neutrino flavor ratio at the source with AR}^ /R}^ taken to be 15% and the 
uncertainties of neutrino mixing parameters given by Eq. ([3]). We have demonstrated that 
the further distinction between z/^ and Vr in such a high-energy range does not improve the 
constraint on the source flavor ratio. This is a consequence of approximate Vf_i — Vr symmetry. 

This work is supported by National Science Council of Taiwan under Grants No. NSC 96- 
2112-M-009-023-MY3 (T.C. Liu and G.L. Lin) and NSC 97-21 12-M-239-002 (M.A. Huang) 
from the National Science Council of Taiwan. 
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